MR Imaging of Reynolds Dilatancy in the Bulk of Smooth Granular Flows 
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Dense granular matter has to expand in order to flow, a phenomenon known as dilatancy. Here 
we perform, by means of Magnetic Resonance Imaging, direct measurements of the evolution of 
the local packing density of a slow and smooth granular shear flow generated in a split-bottomed 
geometry. The dilatancy is found to be surprisingly strong. The dilated zone follow the region 
of large strain rate and slowly spreads as a function of time. This suggests that the local packing 
density is governed by the total amount of local strain experienced since the start of the experiment. 

PACS numbers: 45.70.-n, 46.65.+g, 



Granular media, such as sand, are conglomerates of 
dissipative, athermal particles that interact through re- 
pulsive and frictional contact forces, and that jam into a 
disordered configuration when no external energy is sup- 
plied When slowly sheared, granulates flow, but to 
do so, they must overcome geometrical (steric) hindrance. 
The resulting expansion of the material is referred to as 
Reynolds dilatancy [2|, l3[ . Precise studies of dilatancy are 
hindered by the opaque nature of granular media, and 
by the complexity of slow grain flows, which in experi- 
mental situations never form a simple linear shear profile 
— measurements of average densities therefore 
cannot probe dilatancy. Hence, many basic questions are 
left unanswered. How much do granular material dilate 
in the flowing zone? Does the dilatancy depend on the 
local strain rate? How far does the dilated zone extend? 
Does the packing density evolve over time? 

Here we address these issues by imaging the 3D packing 
density of smooth and slow granular flows by means of 
Magnetic Resonance Imaging (MRI). The flows are gen- 
erated in a split-bottomed shear cell, where the grain flow 
is driven by the rotation of a bottom disc with respect to 
a cylindrical container (Fig. la). In these systems, a wide 
shear zone emanates from the edge of the disc, making 
them well suited for studies of smooth and slow granular 
flows @, 0, H i, El, El El, El, El- In earlier studies, 
the ratio of the filling height H to the radius of the disc 
R s was found to govern the qualitative nature of the flow 
0, @, S, El, E3, El: When H/R s is less than 0.6, the 
central core material rests on and rotates as a solid with 
the disc. The shear zone then reach the free surface, and 
the three dimensional shape of the region of larg e shear 
resembles the cone of a trumpet [1, H, 
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When 

H/R s > 0.75, the flowing zone do not reach the surface 
but forms a "dome" shaped shear zone in the bulk of the 
material [1, El, 11, EH- For intermediate filling depths 
there is a crossover regime [Tl| • 

By imaging the local packing densities in our shear cell 



as a function of time and for a range of filling heights, 
we explore the relation between flow field and dilatancy. 
The simple picture that emerges is that the amount of 
dilatancy grows with the total amount of local strain ex- 
perienced and saturates when the local strain becomes 
of order one. First, we show that the relative change 
in density in the flowing zone is rather strong and satu- 
rates around 10-15 %. Secondly, we find that the dilated 
zone slowly spread throughout the system as time pro- 
gresses. This spread is consistent with the idea that the 
total, local strain experienced since the start of the ex- 
periment governs the amount of dilatancy — to show 
this, we will reconstruct the flow field in the bulk, com- 
bining previously found scaling relations [1, H, H, El] with 
measurements of the velocity at the surface. Finally, and 
in contrast to small filling heights (H/R s < 0.6) where 
the locations of the dilated zone and the shear zone coin- 
cide, for deep filling heights a relatively long-lived tran- 
sient causes the dilated zone to deviate from the late-time 
shear zone. 

Experimental setup — By detecting the signal of pro- 
tons in oil within seeds, MRI has been used to non- 
invasively probe the dynamics of granular materials 
[Ifl,ll3ll6|,ll7|. While various techniques, based on either 
spin-tagging or employing mixtures of grains with dis- 
tinct oil concentration, have been used to track grain mo- 
tion, we are not aware of prior research employing MRI 
to measure packing densities qb. Here, food grade poppy 
seeds sieved to select diameters of one mm were used. 
Multiple horizontal slices were acquired in an interleaved 
fashion to reduce signal saturation effects [18[. Slices 
were 2 mm thick with in-plane resolution of 1.56 mm x 
1.56 mm. Individual seeds are thus not resolved, but the 
MRI signal probes local packing densities. To compen- 
sate for the overall gradients in the signal, we compare 
signal intensity before and after flow has occurred, and 
report the relative change in density Acj) = (j)(t) / (j)(0) — 1. 
We employ the azimuthal symmetry of the system to av- 
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FIG. 1: (Color) (a) Sketch of the split-bottomed shear cell, 
in which a layer of granular material of depth H is driven 
by the rotation of a disc of radius R s with rate (b-d). 
Comparison of the surface flow field uj(r) and the packing 
density in the bulk after N rotations, (b) H/R s = 0.23 and 
N = 2.5; (c) H/R s = 0.41 and N = 4; (d) H/R s = 0.51 
and N = 4. The straight, vertical lines indicate where uj(r) 
reaches 0.1, 0.5 and 0.9, thus indicating the location of the 
shear zone as observed at the free surface. The white curves 
are a theoretical prediction for the center of the shear zone 
based on Eq. [3] [9] . The black curves show examples of the 
edge of the shear zone, based on Eq. [3] and 2) 



erage the data and focus on A0(iV, r, z\ where TV, r and z 
denote the number of rotations of the disc since the start 
of the experiment and the radial and vertical coordinate 
respectively. 

The experimental geometry placed into the MRI ap- 
paratus is a stationary plexiglass drum of inner radius 
95 mm in which the bottom disk of radius R s = 70 mm 
is slowly rotated with a rate Q of approx 1/8 s _1 (see 
Fig. 1). At this rate, the flow velocities are simply pro- 
portional to n 0, 0, S, S, EqI, El EI El El E3 • The cel1 

is filled with poppy seeds and tapped to obtain a well- 
compacted initial condition of filling depth H. A layer of 
seeds is bonded to the side walls and bottom to obtain 
rough boundaries. To study the temporal evolution of 
the density, we have obtained images after N rotations 
of the bottom disc. During the MRI procedure, the disc 
is stopped — starting and stopping the flow has no ap- 
preciable effect on these slow, non-inertial flows. To facil- 
itate direct comparisons to the flow field, we used video 
imaging to obtain the non-dimensional surface velocity 
profile o;(r), defined as the ratio of the time averaged az- 
imuthal velocity to the driving rate Q @, H, El • We focus 
on late times after transients have died out. 

Dilated zone and Shear zone for low filling heights — 
Figure 1 illustrates shear and dilated zones are directly 
related for low filling heights (H/R s < 0.55). First, 
near the surface, the location of the dilated zone and 
the shear zone coincide. Measurements of the surface ve- 
locity profiles uj(r) are in agreement with earlier studies 



@, 13, E3, El El since we find that 

u(r) = 1/2 + 1/2 erf [(r - R c )/W] , (1) 

where the center location R c in consistent with the simple 
scaling law: 



1 - R c /R s = (H/Rs) 



5/2 



(2) 



Figure lb-d shows that the centers of the dilated and 
shear zones coincide at the free and bottom surfaces, and 
that the width of both zones grows similarly with H/R s . 

Second, in the bulk the dilated zone also follow the 
shear zone. The center of the shear zone throughout the 
bulk R c and the vertical coordinate z are related to the 
geometrical parameters H and R s by a scaling argument 
as Q: 



2.5\\0.4 



H - R c (l - Rs/Rdl - (H/R s y- b )) 



(3) 



Fig. lb-d shows that this curve falls into the center of the 
dilated zone. 

Third, data on the width of the shear zone in the bulk, 
obtained by either excavating colored particles [Si I 1 11 or 
MRI and numerical studies of the flow field [10, LL3[ , in- 
dicate that the shear zone rapidly widen with distance 
from the bottom disc as [19[ 
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We have shown examples of such curves in Fig. lb-d, 
and find that they describe the shape of the dilated zone 
well [HI, Overall, our data clearly suggest that for 
H/R s < 0.55 the shear and dilated zones are directly 
related. 

Spread of the dilated zone — Notwithstanding their 
close relation, the shear zone and dilated zone are not 
identical. To clarify the relation between dilatancy and 
flow field further, we show in Fig. 2a the evolution of the 
density profiles obtained along a number of vertical slices 
as a function of the number of rotations N. 

We first focus on the density profiles for fixed N. The 
local flow and strain rate vary substantially throughout 
the shear zone (see below), but in contrast, the dilatancy 
clearly reaches a plateau value throughout the shear zone. 
Packing density and local strain rate are thus not directly 
related. We have found that this maximal dilatancy is 
similar for all our data sets, and does not vary much with 
time and depth, except close to the free surface [20]. The 
amount of dilatancy is substantial, with the densities in 
the plateau 10-15% lower than the non-sheared regions. 
To put this number in perspective, disordered packings 
of spheres span a range in packing densities from 55% 
(Random Loose Packing) to 64 % (Random Close Pack- 
ing), with typical densities of poured grains around 60 
% [2l H2[. Even though poppy seeds are not perfectly 
spherical (see Fig. 2b), so that the precise values of these 
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FIG. 2: (Color) (a) Horizontal profiles of the density changes 
for the data set with H / R s — 0.51 (shown in Fig. Id), as 
functions of N and for five values of z as indicated. N — 
1/2, 1,2,4, ...64, with N —1/2 shown in orange, and N = 64 
in black. The vertical lines indicate estimates of where, after 
N turns, the local strain rate equals one. (b) Snapshot of 
the poppy seeds used in this experiment, (c) Surface veloc- 
ity co(r) (dashed), a fit to the error function form of Eq. ([J) 
(red), and the local strain after one rotation 2-Krd r u (black, 
log scale) . The horizontal lines indicate where the local strain 
7 reaches one after N rotations — colors as in panel (a), (d) 
Illustration of the locations where 7 after N rotations reaches 
one, constructed from the strain at the surface and the scaling 
laws Eq. J3} (dashed line) and Eq. J4). 



numbers may be somewhat different, the observed dila- 
tancy is substantial. 

The most striking feature shown in Fig. 2a is that the 
dilated zone slowly spreads with the number of rotations 
N [23[. As we will now show, the picture that emerges 
is that the total amount of strain, not the strain rate, 
governs the amount of dilatancy — entirely consistent 
with a quasi-static picture of slow grain flows. 

This is illustrated in Fig. 2. Away from the container 
bottom, vertical gradients are small, and we approximate 
the total local strain in the bulk after N turns by 



7 = 2nNrd r 



(5) 



To obtain d r uo, we start from the surface velocity profile 
(black dashed curve in Fig. 2c), to which we fit the error 
function form of Eq. (pQ) (red curve in Fig. 2c) before 
taking the derivative to obtain the local strain after one 
rotation (27rrd r uj) (black curve in Fig. 2c, log scale). The 
horizontal lines in Fig. 2c indicate where the local strain 
at the surface after N rotations reaches one. 

Figure 2d illustrates how we determine the local strain 
field in the bulk from the surface measurements. We 
start with the observation that the functional form of 
uj(r) in the bulk is also an error function [l(J El[]. The 
center location of this shear zone is given by Eq. [3] (dashed 
curve in Fig. 2d). Finally, in the bulk the width of the 



measured width of the shear zone at the free surface we 
can calculate the complete strain rate field. From this 
we can then deduce, for any number of turns iV, where 
in the bulk the local strain reaches one (colored curves 
in Fig. 2d). 

In Fig. 2a, the density profiles for N = 
1/2,1,2,4, ... ,64 are compared to the location in 
the bulk where after N turns 7 = 1. We see that, at 
least qualitatively, the location where the strain reaches 
one coincides with the edge of the dilated zone. Near 
the container bottom, the coincidence is less good, 
but there vertical gradients, deviations from Eq. [4] and 
boundary effects may be substantial. Nevertheless, our 
data provides clear evidence that the spreading of the 
dilated zone is simply caused by the spreading of the 
position where the local strain becomes of order one 
with time. This suggests that there is no steady state 
and that the dilated zone continually expands since the 
strain rate is nowhere strictly zero [231 ]. 

In prior work focussing on the velocity field for similar 
filling hei ghts , no long time evolution was observed @, 
EL EO, Hll[l3|. This does not necessarily disagree with 
the slow dilatancy. To measure velocities accurately, one 
requires grains to move substantially in order to average 
out fluctuations. However, the system is already fully 
dilated after a local strain of order 1, so that transients 
in the flow field caused by changes in the density escape 
detection — essentially, dilatancy happens in the far tails 
of the velocity profiles 

Dilated zone for large filling heights — In Fig. 3a 
we show the surface velocities and density evolution for 
H/R s = 0.77. For such deep filling depths, the shear zone 
emanating from the edge of the disk meet in the mid- 
dle of the cell and form a dome-like structure, and the 
surface remains essentially stationary 0, Eo|, EH, EH • Sim- 
ple scaling laws for the flow field like Eqs. (p]|4]) are not 
known, but the qualitative features of the dilated zone at 
early times also show the dome-like shape 0, Eo|, EH, EH • 
For later times the dilated zone spreads throughout large 
parts of the system, leading to an appreciable elevation 
and corrugation of the free surface. Slow convection, oc- 
curring for large filling heights, leads to the slow growth 



of a dip in the center of the free surface 

In Fig. 3b we show the surface velocities and density 
evolution for H/R s = 0.61. This corresponds to the in- 
termediate regime, where the velocity profile is weakly 
asymmetric [11|, but there is no precession. We recover 
this asymmetry of o;(r), and also note that the center 
of the surface remains stationary during a short tran- 
sient. The bulk density evolution reveals that the 3D 
flow profile is a complex mix between the "trumpet" 
shape observed for small filling heights and the "dome" 
shape observed for large filling heights 



13, El SI- For 
early times, all the dilatancy is concentrated in a dome 
like structure, consistent with the steady surface for early 
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FIG. 3: (Color) Surface velocity profiles and dilatancy maps 
for (a) H/R s = 0.77 and (b) H/R s = 0.61. See text. 



times. For late times, the shear zone evolves to the trum- 
pet shape leading to uo « 1 near the center of the surface. 
We suggest that the symmetry breaking of the surface 
flow profiles and the dome-like dilated zone might be re- 
lated pjj. Also here a dip develops near the free surface 
(the almost vertical shape of this dip for N = 512 is an 
imaging artefact). 

Discussion — Our data uncovers how granular media 
dilate under steady shear [25]: the amount of dilatancy 
grows with the accumulated strain, and saturates when 
this strain becomes of order one. This fully dilated state 
might very well coincide with what in the engineering 
literature is known as the "critical state", referring to 
a final state of grain flow where rheological properties 
stay constant [14j . How general this correspondence is 
remains to be seen, in particular for shear bands near a 
boundary or free surface 0, [l7| . 
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